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Abstract: In the framework of the Constrained MSSM we re-examine the gravitino as 
the lightest superpartner and a candidate for cold dark matter in the Universe. Unlike in 
most of other recent studies, we include both a thermal contribution to its relic population 
from scatterings in the plasma and a non-thermal one from neutralino or stau decays after 
freeze-out. Relative to a previous analysis [1] we update, extend and considerably im- 
prove our treatment of constraints from observed light element abundances on additional 
energy released during BBN in association with late gravitino production. Assuming the 
gravitino mass in the GeV to TeV range, and for natural ranges of other supersym- 
metric parameters, the neutralino region is excluded, except for rather exceptional cases, 
while for smaller values of it becomes allowed again. The gravitino relic abundance 
is consistent with observational constraints on cold dark matter from BBN and CMB in 
some well defined domains of the stau region but, in most cases, only due to a dominant 
contribution of the thermal population. This implies, depending on mg, a large enough 
reheating temperature. If mg > 1 GeV then Tr > 10 7 GeV, if allowed by BBN and other 
constraints but, for light gravitinos, if mg > 100 keV then Tr > 10 3 GeV. On the other 
hand, constraints mostly from BBN imply an upper bound Tr < a few x 10 8 GeV which 
appears inconsistent with thermal leptogenesis. Finally, most of the preferred stau region 
corresponds to the physical vacuum being a false vacuum. The scenario can be partially 
probed at the LHC. 

KEYWORDS: Supersymmetric Effective Theories, Cosmology of Theories beyond the SM, 
Dark Matter, Supersymmetric Standard Model. 
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1. Introduction 

Weakly interacting massive particles (WIMPs) remain the most popular choice for cold 
dark matter (CDM) in the Universe. This is because they are often present in various 
extensions of the Standard Model (SM). For example, neutralinos in softly broken low 
energy SUSY models can be made stable by assuming some additional symmetries (like 
i?-parity in SUSY). Their relic abundance in some regions of the parameter space agrees 
with the value of S^cdm^ 2 ~ 0.1 inferred from observations. This last property is sometimes 
taken as a hint for a deeper link between electroweak physics and the cosmology of the 
early Universe. 

However, extremely weakly interacting massive particles (E- WIMPs), 1 have also been 
known to provide the desired values of the CDM relic density. E- WIMPs are particles 
whose interactions with ordinary matter are strongly suppressed compared to "proper" 
WIMPs, like massive neutrinos and neutralinos, whose interactions are set by the SM 
weak interaction strength a wea k ~ 1CP 38 cm 2 times some factors, like mixing angles for 
the neutralino which are often much smaller than one. For such WIMP their physics is 
effectively entirely determined by the Fermi scale and (in the case of SUSY) by the SUSY 
breaking scale MgusY, which, for the sake of naturalness, is expected not to significantly 
exceed the electroweak scale. 

In contrast, a typical interaction strength of E- WIMPs is suppressed by some large 
mass scale m\, oe-wimp ~ ( m w / m h) 2 &weak- One particularly well-known example is the 
gravitino for which m\ is the (reduced) Planck scale M P = 1/^8ttG n = 2.4 x 10 18 GeV. 
Another well-motivated possibility is axions and/or their fermionic partner axinos for which 
the scale m\ is given by the Peccei-Quinn scale f a ~ 10 11 GeV. Both the axion [3] and 

1 Another name used in the literature is 'superweakly interacting massive particles' [2]. 
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the axino [4] have also been shown to be excellent candidates for CDM. Other possibilities 
involve moduli [5] although they strongly depend on a SUSY breaking mechanism. Thus 
the electroweak scale may have little to do with the DM problem, after all. 

In particular, the gravitino as a supersymmetric partner of the graviton, is present 
in schemes in which gravity is incorporated into supersymmetry (SUSY) via local SUSY, 
or supergravity. As a spin-3/2 fermion, it acquires its mass through the super-Higgs 
mechanism. Since gravitino interactions with ordinary matter are strongly suppressed, it 
was realized early on that the particle was facing various cosmological problems [6, 7]. 
If the gravitino is not the lightest superpartner (LSP), it decays into the LSP rather 
late (~ 10 2-10 sec) and associated electromagnetic (EM) or hadronic (HAD) radiation. 
If too much energy is dumped into the expanding plasma at late times > 1 sec, then 
the associated particles produced during the decay (e.g., a photon in the gravitino decay 
to the neutralino) can cause unacceptable alterations of the successful predictions of the 
abundances of light elements produced during Big Bang nucleosynthesis (BBN) , for which 
there is a good agreement between theory on the one hand and direct observations and 
CMB determinations on the other. Since the number density of gravitinos is directly 
proportional to the reheating temperature Tr, this leads to an upper bound of Tr < 
10 6 ~ 8 GeV [8, 9, 10, 11, 12, 13] (for recent updates see, e.g., [14, 15, 16]). On the other 
hand, when the gravitino is the LSP and stable (the case considered in this work), ordinary 
sparticles will first cascade decay into the lightest ordinary superpartner, which would be 
the next-to-lightest superpartner, (NLSP) which would then decay into the gravitino and 
associated EM and/or HAD radiation. A combination of this and the overclosure argument 
(Q^h 2 < 1) has in this case led to a rough upper bound Tr < 10 9 GeV [10, 17]. 

There are some generic ways through which gravitinos (assumed from now on to be 
the LSP) can be produced in the early Universe. One mechanism has just been described 
above: the NLSP first freezes out and then, at much later times, decays into the gravitino. 
Such a process does not depend on the previous thermal history of the Universe (so long 
as the freeze-out temperature is lower than the reheating temperature Tr after inflation). 
As in our previous work, we will call it a mechanism of non-thermal production (NTP). 
In a class of thermal production (TP) processes gravitinos can also be generated through 
scattering and decay processes of ordinary (s)particles during the thermal expansion of the 
Universe. Once produced, gravitinos will not participate in a reverse process because of 
their exceedingly weak interactions. Analogous mechanisms exist in the axino case [4]. In 
addition, there are other possible ways of populating the Universe with stable relics, e.g. 
via inflaton decay or during preheating [18, 19], or from decays of moduli fields [20]. In 
some of these cases the gravitino production is independent of reheating temperature and 
its abundance may give the measured dark matter abundance with no ensuing limit on Tr. 
In general, such processes are, however, much more model dependent and not necessarily 
efficient [21], and will not be considered here. 

Recently, there has been renewed interest in the gravitino as a stable relic and a 
dominant component of cold dark matter. Thermal production was re-considered in [22, 23] 
while non-thermal production in [2, 24, 25, 26]. In [27] both processes were considered in 
the framework of the Minimal Supersymmetric Standard Model (MSSM) in the context of 
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thermal leptogenesis. In [1] some of us considered a combined impact of both production 
mechanisms in the more predictive framework of the Constrained MSSM (CMSSM) [28]. 
The CMSSM encompasses a class of unified models where at the GUT scale gaugino soft 
masses unify to m 1 / 2 and scalar ones unify to m®. We concentrated on mg in the GeV to 
TeV range, typical of gravity-mediated SUSY breaking, and on TP contributions at large 
T R ~ 10 9 GeV. 

Since all the NLSP particles decay after freeze-out, in NTP the gravitino relic abun- 
dance S7~ TP /i 2 is related to ^nlsp^- 2 ~ the relic abundance that the NLSP would have had 

G 

if it had remained stable - via a simple mass ratio 

^NT P/l 2 = ^^ NLSp/l 2. 
U "1-NLSP 

Note that f2~ TP /i 2 grows with the mass of the gravitino mg. 

The gravitino relic abundance generated in TP can be computed by integrating the 
Boltzmann equation from Tr down to today's temperature. In the case of the gravitino, a 
simple formula for f2~ p /i 2 has been obtained in [22, 23] 

G \W 10 GeV ) \ m d J \lTeV J ' [ ] 

where rrig(fj,) above is the running gluino mass. In [22, 23] it was argued that, for natural 
ranges of the gluino and the gravitino masses, one can have f2~ p /i 2 ~ 0.1 at Tr as high as 
10 9 ~ 10 GeV. 

The problem is that in many unified SUSY models, the number density of stable relics 
undergoing freeze-out is actually often too large. For example, in the CMSSM, the relic 
abundance of the lightest neutralino typically exceeds the allowed range, except in relatively 
narrow regions of the parameter space. This can be easily remedied if the neutralino is 
not the true LSP and can decay further, for example into the gravitino (or the axino [4]). 
Indeed, it is sufficient to take a small enough mass ratio in the formula (1.1) above. Then, 
however, f2~ p /i 2 may become too large because of its inverse dependence on mg, especially 
at high values of Tr, essential for thermal leptogenesis [29, 30], and for mg in the GeV to 
TeV range. 

One may want to suppress the contribution from TP by considering Tr <C 10 9 GeV 
and generate the desired relic density of gravitinos predominantly through NLSP freeze- 
out and decay. This would normally require a larger gravitino mass mg and therefore 
longer decay lifetimes (see below). This, however, can lead to serious problems with BBN, 
as discussed above. Furthermore, late injection of energetic photons into the plasma may 
distort the nearly perfect blackbody shape of the CMB spectrum [31]. 

In a previous paper [1] by some of us, the issue of a combined impact of TP and 
NTP mechanisms of gravitino production, in view of requiring the total gravitino fig/i 2 = 
ri~ TP /i 2 + ri~ p /i 2 ~ 0.1 and of BBN, CMB and other constraints, has been examined in 

(j (j 

the framework of the CMSSM. 

In the CMSSM, the NLSP is typically either the (bino-dominated) neutralino (for 
m\j2 <C mo) or the lighter stau f\ (for m^ 2 S> mo). Assuming natural ranges of 



-3- 



m i/2i m o % a f ew TeV, the whole neutralino NLSP region was found [1] to be excluded by 
constraints from BBN because of unacceptably large showers generated by NLSP decays 
even assuming fairly conservative abundances of light elements. This confirmed the find- 
ings of [27, 25]. On the other hand, the fraction of the stau NLSP region excluded by the 
BBN constraint was found to depend rather sensitively on assumed ranges of abundances 
of light elements. 

However, it was also found in [1] that the stau NLSP region of parameter space where 
Q x h 2 ~ 0.1 was due to NTP alone were in most cases excluded by (mostly) the BBN 
constraints. In other words, for natural ranges of mi/ 2 , a significant component of QqIi 2 
from TP (and thus rather high Tr) must normally be included in studies of gravitino CDM 
in the CMSSM. 

In both the neutralino and the stau NLSP cases, decay products lead mostly to EM 
showers but a non-negligible fraction of them develop HAD showers which can also be very 
dangerous, especially at shorter NLSP lifetimes (< 10 4 sec). Constraints on EM showers 
were analyzed in [14] assuming rather conservative ranges for the abundances of light 
elements but important constraints from HAD showers were not included. A combined 
analysis of constraints on both EM and HAD fluxes was recently performed in [15] with 
much more restrictive (and arguably in some cases perhaps too restrictive) observational 
constraints than [14]. In [1] the EM shower constraint was applied following [14] and 
the HAD one following [15] and thus, out of necessity, using different assumptions about 
allowed ranges of light elements. 

In the present analysis we make a number of improvements. Firstly, we treat both EM 
and HAD showers in a self-consistent way by assuming the same ranges of abundances of 
light elements which we take to be somewhat less restrictive than those adopted in [15]. 
Secondly, the BBN yields are computed with a sophisticated code which simultaneously 
deals with the impact of both EM and HAD showers. In some parts of the parameter space 
this is essential since non-linearities may exist and in such cases EM and HAD activities 
cannot be analyzed separately in computing the abundances of light elements. Thirdly, at 
each point in the CMSSM we compute energy released into all relevant (EM and HAD) 
channels and their hadronic branching ratios (which are typically smaller than the EM ones 
but can vary by a few orders of magnitude) and then use them as inputs into the BBN 
code. Furthermore, in [1] in dealing with EM showers we conservatively only included 
bounds from D/H, Y p ( A He abundance) and 7 Li/H while in constraining HAD showers 
we dropped the lithium constraint. In the current work we include all the above three 
constraints and in addition apply constraints from 3 He/D and e Li/ 7 Li which in some cases 
have the strongest impact on the CMSSM parameter space. All these improvements lead 
to a much more reliable BBN constraint on EM/HAD showers, even though a large part 
of the stau NLSP region still remains allowed. In particular, we improve the upper bound 
of T R < 5 x 10 9 GeV found in [1] to Tr < a few x 10 8 GeV. 

In addition, we have now extended the range of mi/2 to 6 TeV, beyond that considered 
in [1], and found that at very large values of ^ 4 TeV one can find allowed regions 
consistent with Q^h 2 in the interesting range due to NTP alone. We have also found 
relatively confined pockets of the neutralino NLSP parameter space which are consistent 
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with BBN and CMB. We discuss these results below. 

Lastly, in SUSY theories the presence of several scalar fields carrying color and elec- 
tric charge allows for a possible existence of dangerous charge and color breaking (CCB) 
minima [32] - [37] which would render the physical (Fermi) vacuum unstable. Along some 
directions in field space the (tree-level) potential can also become unbounded from below 
(UFB). Avoiding these instabilities leads to constraints on the parameter space among 
which those derived from requiring the absence of UFB directions are by far the most 
restrictive [34]. In the specific cases applicable to the CMSSM, after including one-loop 
corrections, such UFB directions become bounded but develop deep CCB minima at large 
field values away from the Fermi vacuum of our Universe. 

Although the existence of such a dangerous global vacuum cannot be excluded if the 
lifetime of the (metastable) Fermi vacuum is longer than the age of the Universe (in addition 
to having rather unpleasant scatological consequences), such a possibility does place non- 
trivial constraints on inflationary cosmology. One has to explain why and how the Universe 
eventually ended up in the (local) Fermi minimum [37, 39]. 

The effect of the UFB constraints in minimal supergravity models was analyzed in [40] , 
where it was shown that it is the stau region in the CMSSM that is mostly affected. This 
is precisely the region of interest for gravitino CDM in the CMSSM. Consequently, we will 
discuss the impact of these constraints in our analysis. 

As in [1], we will take mg as a free parameter and allow it to vary over a wide range of 
values from O(TeV) down to the sub-MeV range, for which the gravitino (at least those 
produced in TP, see later) would remain cold DM relic. Lighter gravitinos would become 
warm and then (sub-keV) hot DM. We will not address the question of an underlying (if 
any) supergravity model and SUSY breaking mechanism. As in [1], we will mostly focus on 
the 0{ GeV) to 0{ TeV) mass range, as most natural in the CMSSM with gravity-mediated 
SUSY breaking, but will also explore at some level light gravitinos. 

In the following, we will first summarize our procedures for computing fig/i 2 via both 
TP and NTP. Then we will list NLSP decay modes into gravitinos, and discuss constraints 
on the CMSSM parameter space, in particular those from BBN and CMB. Finally, we will 
discuss implications of our results for thermal leptogenesis and for SUSY searches at the 
LHC. 

2. Framework and Procedure 

Unless otherwise stated, we will follow the analysis and notation of the previous paper [1] to 
which we refer the reader for more details. Here we only summarize the main points, while 
below we elaborate on our improved analysis of constraints from BBN and on a previously 
neglected impact of CCB minima and UFB directions. 

Within the framework of the CMSSM we employ two-loop RGEs to compute both di- 
mensionless quantities (gauge and Yukawa coupling) and dimensionful ones (gaugino and 
scalar masses) at the electroweak scale. Mass spectra of the CMSSM are determined in 
terms of the usual five free parameters: the previously mentioned tan/3, m^j and mo, as 
well as the trilinear soft scalar coupling Aq and sgn(p) - the sign of the supersymmetric 
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tanff= 10, /x>0 m G = m tanff-50, fj,>0 m G = 0.2 m 




Figure 1: The plane (roi/2,mo) for tan/3 = 10, mg = too (left window) and tan/3 = 50, = 
0.2mo (right window) and for Ao = 0, /i > 0. The light brown regions labelled "LEP x + " and "LEP 
Higgs" are excluded by unsuccessful chargino and Higgs searches at LEP, respectively. In the right 
window the darker brown region labelled "6 — > S7" is excluded assuming minimal flavor violation. 
The dark grey region below the dashed line is labelled "TACHYONIC" because of some sfermion 
masses becoming tachyonic and is also excluded. In the rest of the grey region (above the dashed 
line) the stau mass bound uif 1 > 87 GeV is violated. In the region "No EWSB" the conditions of 
EWSB are not satisfied. Magenta lines mark contours of the NLSP lifetime tx (in seconds). The 
dotted line is the boundary of neutralino (x) or stau (f) NLSP. 

Higgs/higgsino mass parameter fi. The parameter [i is derived from the condition of elec- 
troweak symmetry breaking and we take \i > 0. We compute the mass spectra with the 
help of the package SUSPECT v. 2.34 [41]. For simplicity we assume i?-parity conser- 
vation even though E-WIMPs, like gravitinos or axinos, can constitute CDM even when 
it is broken. This is because, even in the presence of -R-parity breaking interactions the 
E-WIMP lifetime will normally be very large due to their exceedingly tiny interactions 
with ordinary matter. 

We compute the number density of the NLSP after freeze-out (the neutralino or the 
stau) with high accuracy by numerically solving the Boltzmann equation including all 
(dominant and subdominant) NLSP pair annihilation and coannihilation channels. For a 
given value of mg, we then compute the NTP contribution to the gravitino relic abundance 
via eq. (1.1). In computing the thermal contribution f2~ p /i 2 we employ eq. (1.2). 

After freeze-out from the thermal plasma at t ~ 10~ 12 sec, the NLSPs decay into 
gravitinos at late times which strongly depend on the NLSP composition and mass, on mg 
and on the final states of the NLSP decay. Expressions for Tx = 1/tx> where X denotes 
the decaying particle 2 have been derived in [26, 25]. Given some discrepancies between the 

2 From now on we will denote X — x, fi for brevity. 
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two sources, below (as in [1]) we follow [25]. Roughly, for mg <C mx the lifetime is given 
by 

. /100GeV\ 5 / m d x 2 

rx~10 8 sec £ — . (2.1) 

V m x J V 100 GeV/ V ; 

The exact value of the NLSP lifetime in the CMSSM further depends on a possible relation 
between mg and mi/ 2 and/or mo but in the parameter space allowed by other constraints 
it can vary from > 10 8 sec at smaller mx down to 10 2 sec, or even less, for large m 1 j 2 
and/or mo in the TeV range. 

In fig. 1 we present the gravitino relic abundance and the NLSP lifetime in the usual 
plane spanned by and mo for two representative choices: tan/3 = 10 and mg = mo 
(left window) and tan/3 = 50 and mg = 0.2mo (right window), and for A$ = and fi > 0. 
At small mo and large tan (3 some sfermions become tachyonic, as encircled by a dashed line 
inside the grey region (labelled "TACHYONIC" ) in the right window. Relevant collider and 
theoretical constraints (but not yet those coming from BBN or CMB) are shown. We apply 
the same experimental bounds as in [1]: (i) the lightest chargino mass m^± > 104 GeV, 
(ii) the lightest Higgs mass m h > 114.4 GeV, (iii) BR(B -»• X s j) = (3.34 ± 0.68) x 10~ 4 . 
In addition, now we further impose a stau mass bound m^ > 87 GeV [42] which slightly 
enlarges the grey region beyond the part labelled "TACHYONIC". In this analysis we 
update the top quark mass to the current value of mt = 172.7 GeV [43]. 

To help understanding this and subsequent figures, we remind the reader of some basic 
mass relations. The mass of the gluino is roughly given by m g ~ 2.7mi/2- The mass of 
the lightest neutralino, which in the CMSSM is almost a pure bino, is m x ~ OAmi/2- The 
lightest stau f\ is dominated by tr and well above mz its mass is (neglecting Yukawa 
contributions at large tan/3) roughly given by m 2 ^ ~ m\ + 0.15m^ 2 . This is why at 
mo <C m 1 j 2 the stau is lighter than the neutralino while in the other case the opposite 
is true. The boundary between the two NLSP regions is marked with a roughly diagonal 
dotted line. (In the standard scenario the region of a stable, electrically charged stau relic 
is thought to be ruled out on astrophysical grounds.) Regions corresponding to the lightest 
chargino and Higgs masses below their LEP limits are appropriately marked and excluded. 
Separately marked for tan /3 = 50 is the region inconsistent with the measured branching 
ratio BR(B — ► X s j). (For tan/3 = 10, and generally not too large tan/?, this constraint is 
much weaker and "hides" underneath the above LEP bounds.) In the grey wedge of large 
mo conditions of EWSB cannot be satisfied. Finally, for some combinations of parameters 
the gravitino is not the LSP. We exclude such cases in this analysis. 

Also shown in fig. 1 are contours of the NLSP lifetime tx (in seconds). Their shape 
strongly depends on gravitino mass relation with mo, m\/2 and, because of SUSY mass 
relations, on other parameters which determine mx, but in the cases considered here, at 
small to moderate mo, tx typically decreases with increasing m^ 2 . 

Finally, by comparing fig. 1 with fig. 2 of [1] (where m t = 178 GeV was assumed), we 
can see the sensitivity to the top quark mass. The region disallowed by the LEP Higgs 
mass bound has now broadened from m^ ~ 300 GeV to 400 GeV and at low and 
large mo a wedge inconsistent with EWSB has appeared. There is also a noticeable change 
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is in the pseudoscalar resonance region for tan/3 = 50. The green band moves to smaller 
rrii/2 and higher ttiq for smaller rrit- 

3. Improved BBN Analysis 

NLSP (x or fi) decays after freeze-out can generate highly energetic electromagnetic and 
hadronic fluxes which can significantly alter the abundances of light elements. At longer 
lifetimes tx > 10 4 sec EM constraints are strongest but at earlier times HAD shower 
constraint typically become dominant while the EM one virtually disappears. However, 
even at later times HAD constraints can be important. 

We will evaluate the abundances of light elements produced during BBN in the presence 
of EM/HAD showers and, by comparing them with observations, place bounds on the latter. 
To this end, we need to know the energy ef transferred to each decay channel i = em, had 
and their respective branching fractions to EM/HAD showers Bf as well as the NLSP 
lifetime tx- All the above quantities depend on the NLSP and (with the exception of the 
yield) on its decay modes and the gravitino mass. For the cases of interest (% and f±) these 
have been recently evaluated in detail in [25]. In [1] and below we follow their discussion. 

For the neutralino NLSP the dominant decay mode is % — > Gj. In the CMSSM the 
neutralino is a nearly pure bino, thus % ~ B. The decay % produces mostly EM 

energy. Thus 

2 2 

mi, — m~ 

f x — x G dx ~ i ftii 

t em — 2 m ' em ~ y - 1 ) 

Above their respective kinematic thresholds, the neutralino can also decay via x ~^ 
GZ, Gh, GH, GA for which the decay rates are given in [25, 26]. These processes contribute 
to HAD fluxes because of large hadronic branching ratios of the Z and the Higgs bosons 
(Bfo ad ~ 0.7, B-^ ad ~ 0.9). In this case the transferred energy ef and branching fraction 
Bf each channel are 

ml-m% + ml rfy -> Gk)B£ ^ 

e ^^S^ B U*= rj k = Z,h,H,A, (3.2) 

where 

T tot ~ r ( x ^ G 7 ) + Yl r (x - Gk) . (3.3) 

k 

Below the kinematic threshold for neutralino decays into G and the Z/Higgs boson, 
one needs to include 3-body decays with the off-shell photon or Z decaying into quarks 
for which w \{rn x — tHq) and B^ ad (x — > Gj*/Z* — > Gqq) ~ 10 -3 [25]. This provides a 
lower bound on Bf iad . At larger m x , Higgs boson final states become open and we include 
neutralino decays into them as well. 
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The dominant decay mode of the lighter stau fi is f\ — ► Gr which, as argued in [24, 25], 
contributes basically only to EM showers. Thus 

(3.4) 

2 2 

where the additional factor of 1/2 appears because about half of the energy carried away 
by the tau-lepton is transmitted to final state neutrinos. Our results are not sensitive to 
an order of two variations of this overall pre-factor. 

As shown in [25], for stau NLSP, the leading contribution to HAD showers come from 
3-body decays f\ — > GtZ,Gv t W, or from 4-body decays fi — > Gt^*/Z* — > Grqq. The 
transferred energy ef- and branching fractions i (i = Z, W, qq) are 



and 



B ft _ r(n -» GrZ)^ fl _ r(n -» gM^)g^ fl _ rfo - gigg 7) 

D had, Z — r ' -°/iad, — r ' ^/lad, qq ~ p ' 

J- tot J- tot J- tot 

where 

r tot ~ r (fi Gr) + r(fi GrZ) + r(fi Gu T w). (3.8) 

One typically finds ~ 10~ 5 — 10~ 2 when 3-body decays are allowed and ~ 10~ 6 from 
4-body decays otherwise, thus providing a lower limit on the quantity [25]. (Since the 
process f\ — ► Gv T W in proportional to the fx component of f\, which in the CMSSM is 
suppressed, its contribution to Ttot is likely to be tiny) Given such a large variation in 
Bhadi tne choice (3.6) is probably as good as any other. 

For each point in the parameter space and for a given mg, the partial energies ef 
released into all the channels and their respective branching fractions Bf are passed on 
to the BBN code which computes light element abundances in the presence of additional 
EM/HAD showers. A determination of ef for all the individual hadronic channels is 
necessary since the changes of light element yields are not simply a linear function of ef. 
(Note that in our previous analysis [1] this yield dependence on ef was neglected.) The 
output is then compared with observational constraints. 

We emphasize that our treatment of the branching ratios constitutes a significant 
improvement relative to previous analyses where only a few sample calculations at a fixed 
NLSP mass (e.g., 100 GeV or 1 TeV) and for only a fixed hadronic branching ratio (e.g., 
Bhad = 10~3) were linearly extrapolated to derive BBN yield predictions. This is inaccurate 
for several reasons. Hadronic BBN yields are not simply linear functions of the hadronic 
energy injected, but rather depend in a more complicated way on the energy of the hadronic 
primaries. Furthermore, at early times (r < 10 4 sec) a simple linear extrapolation from one 
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calculation with particular B^ad is impossible due to the interplay between the hadronic 
perturbations and thermal nuclear reactions. Finally, at later times, cancellation effects 
between HAD and EM light element production and destruction processes may occur. All 
these effects may be properly addressed only when for each point in the SUSY parameter 
space an separate BBN calculation is performed. 

At early times r < 10 4 sec limits from BBN on particle decay induced showers come 
from injections of hadrons, i.e., mesons for r < 10 2 sec and nucleons for r > 10 2 sec, with 
EM showers having no effect at such early times. Mesons convert protons to neutrons 
by charge exchange reactions, e.g. tt~ + p —> tt° + n [45], thereby increasing the final 
4 He abundance. Nucleons lead to an increase in the D abundance due to both injected 
neutrons fusing to form D or inducing the spallation of A He and concomitant production 
of D [46]. Injected energetic nucleons may also affect a very efficient 6 Li production for 
r > 10 3 sec [46, 47, 15]. At times r > 10 sec HAD showers are still important in setting 
constraints, unless B^ ad is very small. In addition, EM showers, also lead to distortions of 
the light element abundances by photo disintegrating elements [10]. EM showers typically 
lead to elevated 3 He/D [49, 15] and 6 Li/ 7 Li [50] ratios. (Note that the effects of HAD 
showers have not been considered in ref. [26]). 

In the present analysis we fix the baryon-to-photon ratio r\ at 6.05 x 10~ 10 which 
is consistent with the WMAP result r/ = 6.1±q% x 10" 10 [44]. All processes required for 
an accurate determination of the light element abundances are treated in detail. The 
calculations are based on the code introduced in ref. [47] with the effects of EM showers 
added. Details of this code will be presented elsewhere [48]. (A similar detailed presentation 
can be found in ref. [15].) 

We apply the following observational constraints 



2.2 x 10~ 5 < 


D/H 


< 5.3 x 10 


0.232 < 




< 0.258 


8 x 10- 11 < 


7 Li/H 






3 He/D 


< 1.72 




6 Li/ 7 Li 


< 0.1875. 



Note that in [15] a much less conservative upper bound of 3.66 x 10~ 5 on D/H was 
assumed. On the other hand, the constraint from 6 Li/ 7 Li was not applied. Note also that, 
relative to [1], we now also include important constraints from 3 He/D and e Li/ 7 Li. These 
bounds are typically a factor of ten more stringent than constraints from D alone. Though 
the stellar evolution of 3 He is not well understood, the constraint from the 3 He/D ratio 
is very secure as D is known to always be destroyed in stars whereas 3 He may be either 
destroyed or produced in stars. Thus one may derive an upper limit on the primordial 
3 He/D ratio given by its value in the pre-solar nebula. Less secure is the upper limit 
on the e Li/ 7 Li ratio, as e Li is more fragile than 7 Li. Typical G Li/ 7 Li observations in low- 
metallicity stars fall in the range ~ 0.03 — 0.1 (cf. [51]). 

Finally we mention the constraint from the CMB shape. Late injection of electromag- 
netic energy may distort the frequency dependence of the CMB spectrum from its observed 
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Figure 2: The same as in fig. 1 but with constraints from BBN and CMB superimposed. The 
regions excluded by the various BBN constraint are denoted in violet. The region disallowed by 
D + Y p and additional regions excluded by 3 H e and 6 Li are denoted accordingly by their respective 
names. A solid magenta curve labelled "CMB" delineates the region (on the side of label) incon- 
sistent with the CMB spectrum. In both windows, the dark green bands labelled "T^ = 10 7 GeV" 
and "10 s " denote the total relic abundance of the gravitino from both thermal and non-thermal 
production with denoted reheating temperature is in the favored range, while in the light green 
regions (marked "NTP") the same is the case for the relic abundance from NTP processes alone. 

blackbody shape [9, 31], as recently re-emphasized in [24] and in [1]. In this paper we apply 
our analysis presented in [1] to which refer the reader for all the details. 

In fig. 2 we present in the usual (rrii/2, m o) plane our new and more accurate constraints 
from BBN for the same parameters as in fig. 1. First we note that a robust constraint from 
D/H and Y p excludes basically the whole neutralino NLSP region. 3 This remains true 
so long as > 0(1 GeV) as we will discuss later. This is consistent with the previous 
analysis [1] and also confirms the findings of [27, 25]. Next, the important role played by 
the constraints from 3 He/D and e Li/'Li in excluding additional regions of the (n-i/2) m o) 
is shown explicitly. 

The cases presented in fig. 2 correspond to some of the cases presented in fig. 2 in 
ref. [1] in order to allow a comparison with a previous approximate treatment of the BBN 
constraint. We can see a substantial weakening of the HAD shower bound coming from 
D/H. This is because in ref. [1] a much stronger upper bound on the allowed abundance 
of deuterium was applied, following ref. [15]. Otherwise, the general pattern of excluded 
regions is roughly similar. 

The total gravitino relic abundance consistent with the 2 a range 0.094 < ftcDM^ 2 < 
0.129 [44] (marked Og/i 2 ) is shown in dark green. For comparison, light green regions 

3 However, we do find some limited exception to this, as we discuss below. 
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(marked NTP) correspond to the gravitino relic abundance due to NTP alone in the same 
range. These regions become cosmologically favored when one does not include TP or 
when Tr <C 10 9 GeV. (Note that these regions correspond to ranges of mi/ 2 beyond those 
explored in [26], where only NTP was considered, and were not found there.) In the white 
regions encircled by the dark (light) green bands, the total (NTP-induced) relic abundance 
is too small while on the other side it is too large. Note that the shape of the green bands 
strongly depends on the gravitino mass relation with mo, m-i/2 and/or other parameters [1]. 
As emphasized in [1], at large Tr these two bands of Q^h 2 and S7~ TP /i 2 correspond to very 
different regions of stau NLSP parameter space. 

On the other hand, it is only at such large m^, where the BBN constraint becomes 
much weaker due to a much shorter lifetime, that we find some cases where NTP alone can 
be efficient enough to become consistent with preferred range of CDM abundance. This 
can be seen in the left window of fig. 2. 

We also note that the constraint from not distorting the CMB spectrum (magenta line) 
seems generally less important than that due to BBN [52]. 

It is thus clear that, so long as Tr < a fewlO 8 GeV, one finds sizable regions of rather 
large m\i2 and much smaller rriQ consistent with the preferred range of CDM abundance. 
Unless one allows for very large m-i/2 ^4TeV, a substantial (and, in fact, dominant) TP 
contribution to f^g/i 2 are required. 

4. False Vacuua 

A complete analysis of all the potentially dangerous CCB and UFB directions in the field 
space of the CMSSM, including the radiative corrections to the scalar potential in a proper 
way, was carried out in ref. [34]. As we commented in the Introduction, the most restrictive 
bounds are from the UFB directions, and therefore we will concentrate on them below. 

In the CMSSM there are three UFB directions, labelled in [34] as UFB-1, UFB-2 and 
UFB-3. It is worth mentioning here that the unboundedness is only true at tree level since 
radiative corrections eventually raise the potential for large enough values of the fields. 
Still these minima can be deeper than the usual Fermi vacuum and thus dangerous. The 
UFB-3 direction involves the scalar fields {H u , ui i , . , } with i ^ j and thus leads also 
to electric charge breaking. Since it yields the strongest bound among all the UFB and 
CCB constraints, and for future convenience, let us briefly give the explicit form of this 
constraint. 

By simple analytical minimization of the relevant terms of the scalar potential it is 
possible to write the value of all the z^, , fields in terms of H u . Then, for any value 
of \H U \ < Mgut satisfying 




9 ,2 +9 2 2 2A e . 
the potential along the UFB-3 direction is simply given by 



(4.1) 



M ,„2 , 2 , 2 m„ , 2771 



Li 



U UF B-3 = (m 2 Hu +mi t )\H u \ z + f^{m z L . + m 2 +m z L J\H u \- . (4.2) 

A ej 9 +92 
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Figure 3: The same as fig. 2 but with UFB constraints (solid blue line and UFB label plus a big 
arrow) added. For < 5 TeV and small mo the UFB constraints disfavor the stau NLSP region 
that has remained allowed after applying the BBN and CMB constraints. 

Otherwise 




In eqs. (4.2) and (4.3) X ej denotes the leptonic Yukawa coupling of the jth generation. 
Then, the UFB-3 condition reads 

Vufb-3(Q = Q)> V Fcrmi , (4.4) 

where V~Fermi = — I (g 12 + (y 2 — v 2 ^ 2 , with v u ^ = (H® d ), is the Fermi minimum eval- 
uated at the typical scale of SUSY masses. (Normally, a good choice for Msusy is 
a geometric average of the stop masses.) The minimization scale Q is given by Q ~ 
max( Xtop\H u |, Msusy)- With these choices for Q and Msusy the effect of the one-loop 
corrections to the scalar potential is minimized. Notice from eqs. (4.2) and (4.3) that 
the negative contribution to Uufb-3 is essentially given by the vr? Hu term, which in many 
cases can be large. On the other hand, the positive contribution is dominated by the term 
oc 1/X ej , thus the larger A e the more restrictive the constraint becomes. Consequently, 
the optimum choice of the e-type slepton is the third generation one, i.e. ej = f. 

Moreover, since the positive contribution to Uufb-3 is proportional to m|, the potential 
will be deeper in those regions of the parameter space where the staus are light and the 
condition (4.4) is more likely to be violated. For this reason, the cases with stau NUSP are 
typically more affected by the UFB constraints [40] . 

In fig. 3 we present the cases displayed above in fig. 2 but now in addition we mark 
the regions corresponding to a one-loop corrected UFB-3 direction becoming the global 
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Figure 4: The same as fig. 3 but for fixed gravitino mass, = lOGeV and tan (3 = 10 (left 
window) and = 100 GeV and tan/3 = 50 (right window). 

CCB minimum. These are encircled by a solid blue line and marked "UFB" and a big 
arrow. We can see that, unless is excessively large (mi/2 > 4TeV for tan/3 = 10), in 
both cases the whole previously allowed (white) and cosmologically favored (green) regions 
correspond to a false vacuum. As grows with increasing m-i/2, the UFB constraint 
becomes weaker and eventually disappears. 

As stated in the Introduction, one cannot exclude the possibility that the color and 
electric charge neutral (Fermi) vacuum that the Universe ended up in after inflation is not 
a global one but merely a long-lived local minimum. As discussed in [38, 39], this however 
often puts a significant constraint on models of cosmic inflation. The point is that, at the 
end of inflation, the Universe was very likely to end up in the domain of attraction of the 
global minimum which, even at high temperatures (Tr ~ 10 7-9 GeV), could well have been 
a CCB one. Preventing such situations leads also to constraints on the SUSY parameter 
space. 

For this reason, while the UFB regions presented in fig. 3 (which correspond to the 
Fermi vacuum being a local minimum) cannot be firmly excluded, should SUSY searches at 
the LHC find sparticles with masses indicating such a false vacuum, valuable information 
may be gained about the state of the Universe and about early Universe cosmology. 

In fig. 4 we present two cases with a fixed mg. In the left window tan/3 = 10 and 
rag = 10 GeV while in the right one tan/3 = 50 and mg = 100 GeV. (The "cross" case 
of tan/3 = 10 and = 100 GeV is excluded by a combination of collider and BBN 
constraints.) While now BBN constraints have become somewhat weaker due to smaller 
wig (and therefore smaller tx), the whole neutralino NLSP region is again ruled out as 
well as part of the stau NLSP region corresponding to smaller (and therefore smaller 
mfj and hence larger tx)- Most of the remaining stau NLSP region in both windows 
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corresponds to the Fermi vacuum being a local minimum. 

As mentioned earlier, for one case (tan (3 = 50 and = 10 GeV, not shown here) 
we have found a confined pocket in the neutralino region around m^ — 2.5 TeV and 
mo ~ 5 TeV, close to the region of no EWSB, which is allowed by our BBN and CMB 
constraints. However, we consider it to be an odd exception to the rule rather than a 
typical case. 

Finally, in figs. 3 and 4 for T R = 10 8 GeV one finds m l j 2 < 2 TeV in order to remain 
consistent with f2g/i 2 in the observed range and with BBN constraints. Gaugino mass 
unification relations then imply rrig < 5.4 TeV. With increasing Tr this upper bound goes 
down but we still expect it to be considerably higher than the upper limit m g < 1.8 TeV 
claimed in ref. [27] for Tr = 3 x 10 9 GeV. The difference may be due to the considerably 
less conservative assumptions about the primordial abundances of light elements in ref. [27] 
which followed ref. [15]. 

5. Tr versus mg 

We now extend our analysis to smaller mg down to less than 1 MeV. In general thermal 
relics with masses as low as some lOkeV can constitute cold DM. We note that such small 
values are rather unlikely to arise within the CMSSM with the gravity-mediated SUSY 
breaking scheme where mg is expected to lie in the range of several GeV or a few TeV, 
as mentioned earlier. In other SUSY breaking scenarios mg can often be either much 
smaller or much larger than this most natural range. For example, in models with gauge- 
mediated SUSY breaking [53] the gravitino can be extremely light mg = O(eV). On the 
other hand, in models with anomaly-mediated SUSY breaking [54] gaugino masses are 
typically of order 10 to 100 TeV. These comments notwithstanding, in a phenomenological 
analysis like this one, we therefore think it is still instructive to display more explicitly 
the dependence between mg and cosmological constraints from BBN and CMB and the 
ensuing implications for the maximum Tr. 

On the other hand, in the following we will not apply the UFB constraint. This is 
not only motivated by the reasons discussed above but, more importantly, because the 
constraint strongly depends on the full scalar potential which in turn depends on the field 
content of the model. So far we have assumed minimal supergravity but in other SUSY 
breaking scenarios several new scalars are present which are likely to lead to very different 
UFB constraints. 

In the left windows of figs. 5 and 6 we plot the total gravitino relic abundance f2g/i 2 
(solid line), its TP contribution tt% p h 2 (dot-dashed lines) and its NTP part n% TP h 2 (dotted 
line) as a function of mg, for tan/3 = 10, Aq = 0, [i > and m 1 / 2 = 500 GeV and for 
several choices of Tr. In fig. 5 we take mo = 200 GeV (x NLSP) and in fig. 6 mo = 50 GeV 
(fi NLSP). In both cases the contribution f2~ TP /i 2 from NTP provides a lower limit to the 

G 

total gravitino relic abundance ^g/i 2 while ^g P ^ 2 varies with Tr. 

In the right windows of figs. 5 and 6 we plot the maximum value of Tr consistent 
with 0.094 < Slg/i 2 < 0.129 versus m^, for the same choices of other parameters as in 
the respective left windows. We can see how the strong constraints from BBN and then 
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Figure 5: Left window: The total gravitino relic abundance fig/i 2 (solid lines) as a function 
of the gravitino mass mg for tan/3 = 10, Aq = 0, /i > and for the point mi« = 500 GeV, 
too = 200 GeV (x NLSP). Thermal production contribution (dot-dashed lines) to figft 2 is shown 
for different choices of the reheating temperature (Tr = 10 9 , 10 7 , 10 5 GeV), while the non-thermal 
production one (dotted line) is marked by NTP. The horizontal green band shows the preferred 
range for SIcdm^ 2 (marked WMAP). Right window: The highest reheating temperature (blue line) 
versus mg such that the relic density constraint is satisfied for the same choice of parameters as in 
the left window. The colored regions are excluded by BBN (violet), CMB (right side of magenta 
line), and the gravitino not being the LSP. We can see that the sub-GeV gravitino, Tr as small as 
10 5 GeV are sufficient to provide the expected amount of DM in the Universe. 



CMB only affect larger in the GeV range or more. Sub-GeV gravitino mass leaves 
the CMSSM almost unconstrained by the above constraints. In particular, the neutralino 
NLSP region becomes for the most part allowed again. Increasing mg reduces the effect 
of TP. This is because it becomes harder to produce them in inelastic scatterings in the 
plasma. On the other hand, at some point the bounds from BBN and CMB eventually put 
an upper bound on Tr. We examined a number of cases, including various gravitino mass 
values but could not find consistent solutions above an upper limit of 

T R < a few x 10 8 GeV. (5.1) 

No values of Tr exceeding the above values were also found by considering Aq = 
±1 TeV, in addition to our default value of Aq = 0. When = 1 TeV, the regions excluded 
by constraints from Higgs mass bound and due to a tachyonic region become larger. In 
particular, for tan (5 = 50, the Higgs mass constraint extends to 700 GeV and the tachyonic 
region increases to some mo = 400 GeV and mi/2 = 1000 GeV, while the constraint due 
to B — > X s j becomes weaker and is burried under the Higgs mass constraint. When 
Aq = — ITeV, the Higgs mass constraint become weaker but the tachyonic region become 
even larger and extends to m^ 2 = 1100 GeV for tan/? = 50. For both values of Aq, the 
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Figure 6: The same as fig. 5 but for the different point m 1 / 2 = 500 GeV and mo = 50 GeV (?i 
NLSP). 



neutralino NLSP is still diallowed, while in stau NLSP region some allowed regions remain, 
similarly to the case of Aq = 0. 

By comparing the left and the right windows of figs. 5 and 6, we can again see that the 
NTP contribution alone is normally not sufficient to provide the expected range 0.094 < 
fl^h 2 < 0.129. On the other hand, assuming as small as 100 keV, for the gravitino to 
provide most of cold DM in the Universe, implies that Tr as small as 10 3 GeV is sufficient 
for TP to provide the expected amount of cold DM in the Universe. 

Finally, we comment on the interesting possibility that gravitino relics may not be all 
cold but that a fraction of them may have been warm at the time of decoupling. This is 
because in the case of relics like gravitinos produced in thermal processes at high Tr, their 
momenta exhibit a thermal phase-space distribution while gravitinos from NLSP freeze- 
out and decay have a non-thermal distribution. As a result, even though both populations 
are initially relativistic, they red-shift and become non-relativistic at different times and 
may have a different impact on early growth of large structures, CMB and other observable 
properties of the Universe. This "dual nature" of gravitinos was investigated early on in the 
framework of gauge-mediated SUSY breaking scenario (in which gravitinos are light, in the 
keV range) [55] where the thermal population was warm while the non-thermal one was 
providing a "volatile" component characterized by a high rms velocity i> rms . The scheme 
was originally explored in [56] in the case of light axinos. Depending on the axino mass and 
other properties, they can provide a dual warm-hot distribution [57] or warm-cold one [4]. 
More recently, it was found [58] that late charged (stau) NLSP decays could suppress the 
DM power spectrum if NLSP decays contributed some O(10%) of the total DM abundance 
and tx ~ 10 7 sec. Comparing with our results we conclude that the effect is probably 
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marginal in the CMSSM. A similar effect was also studied in the case of supergravity [59] 
and in a more general setting in neutral heavy particle decays [60]. Large t> rms may lead 
to the damping of linear power spectrum and reducing the density of cuspy substructure 
and concentration of halos, which have been considered to be potential problems for the 
standard CDM scenario. A contribution of such a warm (or hot) component of dark 
matter may, however, be strongly limited by current and future constraints from early 
reionization [61]. In the model studied here gravitinos from TP would provide a cold 
component while those from NTP could be a warm one. A more detailed investigation 
would be required to assess constraints on, and implications of, this mixed warm-cold relic 
gravitino population in the presented framework. 

6. Summary 

We have re-examined the gravitino as cold dark matter in the Universe in the framework of 
the CMSSM. In contrast to other studies, we have included both their thermal population 
from scatterings in an expanding plasma at high temperatures and a non-thermal one from 
NLSP freeze-out and decay. In addition to the usual collider constraints, we have applied 
bounds from the shape of the CMB spectrum and, more importantly, from light elements 
produced during BBN. The implementation of the last constraint has in the present study 
been considerably improved and also updated ranges of light element abundances have 
been used but basically confirm and strengthen our previous conclusions [1]. The neutralino 
NLSP region is not viable, while in large parts of the stau NLSP domain the total gravitino 
relic abundance is consistent with the currently favored range. Unless one allows for very 
large mi/2 ^ 5TeV, for this to happen a substantial contribution from TP is required 
which implies a lower limit on Tr. For example, assuming heavy enough gravitinos (as 
in the gravity-mediated SUSY breaking scheme), mg > 1 GeV leads to Tr > 10 7 GeV (if 
allowed by BBN and other constraints). In a more generic case, if mg > 100 keV then 
T R > 10 3 GeV. Generally, for light gravitinos (mg < 1 GeV) BBN and CMB constraints 
become irrelevant because of NLSP decays taking place much earlier. On the other hand, 
the above constraints imply an upper bound (5.1), which appears too low for thermal 
leptogenesis, as already concluded in [1]. 

Finally, we have shown that in most of the stau NLSP region consistent with BBN 
and CMB constraints the usual Fermi vacuum is not the global minimum of the model. 
Instead, the true vacuum, while located far away from the Fermi vacuum is color and 
charge breaking. 

Implications for SUSY searches at the LHC are striking. The standard missing energy 
and missing momentum signature of a stable neutralino LSP is not allowed in this model, 
unless mg < 1 GeV. Instead, the characteristic signature would be a detection of a massive, 
(meta-)stable and electrically charged particle (the stau). It is worth remembering that 
such a measurement would not be a smoking gun for the gravitino dark matter since in 
the case of the axino as CDM the stau NLSP (as well as neutralino NLSP) is typically 
allowed as well [4]. Finally in some cases it may be possible to accumulate enough staus to 
be able to observe their decays into gravitinos [62] and/or to distinguish them from decays 
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into axinos [63]. SUSY searches at the LHC open up a realistic possibility of pointing 
at non-standard CDM candidates and additionally revealing the vacuum structure of the 
Universe. 
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Erratum 

Following a recent paper of Pradler and Steffen [64] a formula for the thermal pro- 
duction of gravitino in ref. [23] (Bolz, et al.) has been corrected. In addition, we have 
corrected a numerical error in our routine computing a s at high temperatures. As a conse- 
quence, the regions of ^g/i 2 (green bands) in all the figures have shifted to the left, towards 
smaller m 1 / 2 relative to the previous version, which in turn has led to improving the upper 
bound on the reheating temperature by about an order of magnitude. The new bound is 
T R < a few x 10 s GeV. 

We thank J. Pradler and F. Steffen for checking our results and informing us about a 
discrepancy with theirs. 
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